Decay rates and survival probabilities in open quantum systems 
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We provide the first statistical analysis of the decay rates 
of strongly driven 3D atomic Rydberg states. The distribu- 
tion of the rates exhibits universal features due to Anderson 
localization, while universality of the time dependent decay 
requires particular initial conditions. 

PACS numbers: 05.45. Mt, 32.80Rm, 42.50Hz, 72.15Rn 

The macroscopic transport properties of classical and 
quantum systems sensitively depend on the dynamics of 
their microscopic constituents. Whereas regular Hamil- 
tonian dynamics implies the quasiperiodic confinement of 
phase space density, chaotic motion and disorder generi- 
cally induce diffusion-like probability transport exploring 
all accessible phase space at sufficiently long times. In 
low-dimensional, disordered quantum systems Anderson 
localization efficiently inhibits diffusive transport through 
quantum interference, as first predicted by Anderson for 
the charge transfer across disordered solid state samples. 
Later on it was realized that chaos can substitute for 
disorder in Anderson's scenario, leading to dynamical lo- 
calization [1,2] of quantum probability transport. How- 
ever, as a further complication, quantum systems with 
classically non-integrable Hamiltonian dynamics (or, in 
the absence of a clear classical analog, with strongly cou- 
pled degrees of freedom) most often do not feature purely 
chaotic but rather mixed regular chaotic dynamics: be- 
sides the chaotic component of phase space the latter also 
comprises regular regions separated from the chaotic do- 
main by fractal structures like cantori [3]. Hence, dy- 
namically localized quantum transport may be amended 
by quantum tunneling from regular regions or by semi- 
classical localization in the vicinity of partial phase space 
barriers [4]. 

A rather robust measure of transport is the transmis- 
sion or decay probability of some initial probability distri- 
bution across or from a confined region of phase space. 
Whereas the transmission problem immediately suggests 
a scattering approach, decay is most conveniently de- 
scribed by the norm || ■ U \\ of the time evolved 
initial state |^o) upon projection Va on the phase space 
domain we want to characterize. If is confined by an 
absorbing boundary, attached to some leads which allow 
for transport to infinity, or coupled to some continuum 
of states, then || Vn ■ f/|(/>o) II will decrease monoton- 
ically in time, for generic |(/)o) [5-8]. Given the spec- 
tral decomposition of the Hamiltonian, this finite leak- 
age from ri can be accounted for by associating non- 



vanishing (exponential) decay rates Vj with eigenstates 
\ej) with non- vanishing support on fJ. Given some ini- 
tial \4>o) = Vn I'/'o), the survival probability within the 
domain Q, boils down to [6] 

/'surv(i) =11 ■ U \M \h J2 ""3 exp(-r,i) . (1) 

j 

All information on the decay process is now encoded [9] 
(i) in the set of decay rates - which is a global spectral 
property of the problem - and (ii) in the expansion co- 
efficients Wj — |(0o|e)P of the initial wave packet in the 
eigenbasis, what is a specific representation of the initial 
condition, and a local property to the extent that \(f>o) 
is localized. With these premises, our initial assertion 
on the crucial importance of microscopic dynamics for 
macroscopic transport leads to the guiding question of 
our present contribution [9,10]: How do classically mixed 
regular-dynamics affect these global and local spectral 
properties, and, consequently, the time dependence of 
^surv? Whilst, recently, a considerable corpus of liter- 
ature has addressed this problem within the context of 
simple models of mesoscopic transport [7,10,11], we shall 
here underpin our rather general answer by a numerically 
exact treatment of highly excited three dimensional Ry- 
dberg states under strong microwave driving [12]. These 
are paradigmatic objects to (theoretically and experi- 
mentally) study [6,9,12,13] quantum probabihty decay in 
the presence of tunneling as well as semiclassical and dy- 
namical localization, and allow a natural and clean defi- 
nition of Psurv and the associated projection Vn, without 
any approximations. 

Let us start with a brief reminder of dynamical local- 
ization in periodically driven hydrogen Rydberg atoms, 
to establish the analogy between the general transport 
problem depicted above and our specific choice of its 
physical realization: The atoms are initially prepared in 
an unperturbed atomic eigenstate \4>o) = ["-o "mo) la- 
beled by the principal quantum number no, the angular 
momentum Iq, and the angular momentum projection 
niQ. For linearly polarized coherent driving (which we 
treat in dipole approximation) of amplitude F at fre- 
quency uj, niQ is a good quantum number inherited from 
the rotational symmetry of the problem, and the atomic 
excitation process proceeds through the subsequent ab- 
sorption and/or emission of photons from/into the field, 
along a sequence of near resonantly coupled states \n £). 
Provided F is large enough compared to the attractive 
force exerted on the electron by the nucleus, the clas- 
sical dynamics of the driven Rydberg electron initially 
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launched along a Kepler orbit turns chaotic, leading to 
an efficient energy absorption also on the quantum level, 
which finally leads to the ionization of the atom. Hence, 
we define the electronic survival probability as the re- 
maining bound state population after an interaction time 
t, with Vn the projector on the bound state component 
of Hilbert space. The Tj in Eq. (1) then become the ion- 
ization rates of the Floquet eigenstates \ej) of the atom 
in the field, and the LUj the associated overlaps (averaged 
over one field cycle [6]) of |(/>o) with \€j). With some 
theoretical and numerical effort, all these quantities can 
be extracted from the fundamental Floquet eigenvalue 
problem, for typical parameter values of state of the art 
laboratory experiments [13] . 
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FIG. 1. Distribution of the ionization rates F of mi- 
crowave-driven (a) ID and (b) 3D mo = Rydberg states 
of atomic hydrogen, for throe different values of the localiza- 
tion parameter £ — 0.25 (diamonds), 0.5 (stars), 1.0 (circles), 
and 2.0 (pyramids; only in (a), because of the high cost of 
3D calculations at the corresponding high field amplitudes). 
The distributions were generated, at fixed C, by sampling 
the Floquet spectra over the frequency (and corresponding 
amplitude, F) ranges u/2-k = 13.16 ... 16.45 GHz (ID) and 
35.5 ... 36.1 GHz (3D), within a Floquet zone of width uj cen- 
tered around no = 70 (3D) and 100 (ID), respectively. 

The atomic ionization process is mapped on the An- 
derson model through the atomic localization parameter 
£ = ^/r] [2], with ^ the localization length in units of the 
photon energy, rj is the ionization potential of \(t>o), mea- 
sured in units of lu, and is completely determined by no, 
due to the angular momentum degeneracy of the hydro- 
gen spectrum. Provided that no through its correspon- 
dence with the principal action of the classical dynamics 
- is chosen within the chaotic component of phase space 
(possibly garnished by nonlinear resonance islands and 
remnants of classically regular motion immcrgcd in the 
"chaotic sea"), rj can be conceived as a measure of the 
extension of the domain of complex transport along the 
energy axis. If we identify rj with the sample length in 
the original Anderson problem and average over different 



realizations of a fixed value of £ <C 1, quantum suppres- 
sion of transport through exponentially localized eigen- 
functions over the quasi-resonantly coupled states on the 
energy axis [2,7,14] imphes a distribution P(r) ~ l/F of 
the ionization rates entering Eq. (1), in accord with the 
current theory [11] on open Anderson insulators. 

Statistically independent realizations of a fixed local- 
ization parameter can be generated by simultaneous vari- 
ation of F and uj, due to the relation jC ~ 6.66F^no/cj^/^ 
from the original theory on dynamical localization [2] 
in periodically driven, one dimensional hydrogen atoms. 
Whilst this theory has no quantitative predictive power 
[6,12], it is nonetheless known to give qualitatively cor- 
rect predictions, even for 3D Rydberg states with not too 
large values of the angular momentum quantum number 
£. In the sequel we therefore use this quantity such as to 
guide our choice of F and lu values. 

Fig. 1 shows the distribution of ionization rates F of 
real 3D hydrogen Rydberg states and of a one dimen- 
sional model atom (where the Rydberg electron is con- 
fined to ID configuration space z > defined by the 
polarization axis, with the Coulomb singularity at the 
origin [6,9]) exposed to a microwave field, within a Flo- 
quet zone (of width uj on the energy axis [6,13]) centered 
around the ng = 70 (3D) and ng = 100 (ID) Rydberg 
manifold. The different values of £ = 0.25... 2.0 are 
realized by sampling the spectra over frequency ranges 
a;/27r = 35.5 . . . 36.1 GHz (3D), and 13.16 . . . 16.45 GHz 
(ID), respectively, and adjusting F correspondingly. Due 
to the dramatically enhanced spectral density of the 3D 
as compared to the ID problem (a consequence of the ad- 
ditional degree of freedom labeled by i) only ten equidis- 
tant w-values are needed to generate the 3D data (with 
a total sample size of approx. 25000 eigenstates), in con- 
trast to 500 in the ID case (up to 100000 eigenvalues). 

For small localization parameter C = 0.25, both, the 
ID model atom as well as the real object do exhibit ion- 
ization rates distributed according to the F"^ law, over 
approx. seven orders of magnitude from F ~ 10~^^ to 
F ~ 10~^ a.u., in nice agreement with the prediction for 
decay from a disordered solid [11]. However, as the local- 
ization parameter is increased (by systematically increas- 
ing F over the entire frequency range indicated above) 
to C = 1.0, we observe a depletion of the distribution at 
small rates, balanced by a decrease of the slope of P{T) 
in an intermediate range F ~ 10^^" . . . 10^^ a.u. Such 
behavior is clearly incompatible with the simple picture 
of exponentially localized probability densities along the 
energy axis (due to dynamical localization of the eigen- 
states over the chaotic phase space component) [2,7,14], 
and suggests the presence of alternative transport mech- 
anisms which compete with dynamically localized diffu- 
sion [9,10]. As already anticipated in the introduction, 
such observation does not come as a surprise, since pe- 
riodically driven Rydberg states exhibit a mixed regu- 
lar chaotic phase space. Indeed, it is known from ex- 
perimental as well as exact theoretical/numerical studies 
[9,12,15] of strongly perturbed atomic Rydberg systems, 
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and, more recently, from mesoscopic systems with de- 
cay [10], that the oigcnstates of such quantum systems 
can be classified according to their localization [9,10,16] 
(i) on regular and/or elliptic regions, (ii) along remnants 
of regular motion immerged in the chaotic sea ( "separa- 
trix/hicrarchical states" living on "can-tori"), and (iii) in 
the chaotic domain of phase space. Furthermore, elliptic 
regions as well as remnants of invariant tori are rather 
robust under changes of some control parameter such as 
F in our present case [9] . 
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FIG. 2. (a) Expansion coefficients Wj of the 3D atomic ini- 
tial state jno = 70, io = 0, mo = 0} in the Floquet eigenstates 
\ej), for £ = 1.0 and ij/2-K = 35.6 GHz, vs. the associated 
decay rates Tj . (Qualitatively similar results are obtained for 
the one dimensional model atom, as well as for microwave 
driven alkali Rydberg states [17].) (b) Same as (a), but for 
initial states prepared at the sites 900 (plusses), 999 (crosses), 
1000 (diamonds clustered along the diagonal) of a ID Ander- 
son model of sample length A'^ = 1000 (uniformly distributed 
on-site potentials Ki; with jKiJ < 1.5, constant coupling 
Vn^,ni±i = 1, and an absorbing boundary at one end, in- 
troduced by adding an imaginary part —ij = — t X 0.31 to the 
on-site potential at rii — 1000; results for 10 realizations of 
the potential are gathered together) . Only in the special case 
of the initial state placed right at the edge ni = 1000 of the 
sample is there a one-to-one relation between decay rates and 
expansion coefficients! 

Accordingly, the localization properties as well as the 
decay rates of the associated eigenstates tend to be less 
sensitive when changing F, whereas those states localized 
in the chaotic domain, with their decay rates essentially 
determined by C, exhibit rapidly increasing rates as F 
is increased. Hence, we attribute the deviation from the 
law at large values of F in Fig. 1 to the rapid increase 
of the decay rates of class (iii) states with C (correspond- 
ing to the effective destruction of dynamical localization), 
whereas class (i) and (ii) states accumulate at intermedi- 
ate F-values at the largest £-values considered here (we 
verified this assertion by inspection of the Husimi phase 



space representations of the Floquet eigenstates of the 
ID problem [9]). At small values of £, the latter exhibit 
extremely small ionization rates below F ~ 10"""^^ a.u. 
(or even below F ~ 10~^^ a.u., this is below our nu- 
merical precision, hence the plot does not extend into 
that domain), what leads to isolated, extremely narrow 
resonances in experimental or numerical photoionization 
cross sections [9,15]. Therefore, in this parameter regime, 
P(F) is dominated by the decay rates of the dynamically 
localized class (iii) states, in agreement with the Ander- 
son scenario. Let us also comment on the apparently 
distinct weight of class (i) and (ii) states in the ID and 
the 3D situation, which is born out in Fig. 1: the de- 
pletion of the low-F part is less pronounced for the real 
atom, and equally so the region of decreased slope at in- 
termediate F-values. This is due to the contribution of 
high £ states (at fixed mo) to the rate distribution - un- 
available in ID - which exhibit vanishing decay rates and 
do not contribute to the plotted distribution at small C 
As C is increased, these states contribute to P(F) in the 
range F < 10~^^ a.u., such as to compensate for the re- 
shuffling of low-£ states from small to intermediate F's. 
This smooths the 3D distribution as compared to the ID 
result. 

Hence, we have seen how chaotic transport can mimic 
signatures of disorder in the quantum decay rate distri- 
bution, and induce amendments thereof due to peculiar 
structures of mixed regular chaotic dynamics. As obvi- 
ous from Eq. (1), P(F) does not determine the (exper- 
imentally often very easily accessible [6,9,13]) survival 
probability alone, but only in conjunction with the ex- 
pansion coefficients Wj. Starting from a general form 
P(F) ~ r~", many authors infer a "universal" time de- 
pendence Psurv ^ 1/^^"", leaning on the additional as- 
sumption Wj ~ Tj [7,8]. The latter is nonetheless only 
justified if the initial state \4>o), which represents a local 
probe of the spectrum, is localized close to the boundary 
of f2 - which is the "edge of the sample" in a language 
adapted to the Anderson scenario, and the ionization 
threshold n — > oo in the atomic problem. For a generic 
choice of {(po) "within the sample" , i.e. somewhere within 
ri, this basic assumption is utterly misleading, as illus- 
trated in Fig. 2. In Fig. 2 a), we show the distribution 
of the Wj and of the corresponding F^, on a double log- 
arithmic plot, for ]^o) = I'T^o = 70 £q = 0) in driven 
3D hydrogen. The essentially uncorrelated Wj and Fj 
nicely correlate with the earlier finding that phase space 
localization properties and decay rates of strongly driven 
Rydberg atoms are not unambiguously related [9], and 
clearly devalidate the above proportionality assumption. 
On the other hand, it can easily be shown within a simple, 
one dimensional Anderson model how the generic situa- 
tion of Fig. 2 a) continuously evolves into the situation of 
strongly correlated Wj and Fj, as illustrated in Fig. 2 b): 
There, we calculated the decay rates of the eigenstates 
of a ID tight binding Hamiltonian of length N = 1000, 
with diagonal disorder and constant off-diagonal terms. 
Decay is introduced in this model by absorbing bound- 
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ary conditions at one end. whereas a reflecting wall is 
imposed at the other. If \4>q) is placed well inside the 
sample, expansion coefficients and decay rates are uncor- 
related, though the distribution collapses onto a straight 



line Wi 



Tj as 



t>o) is shifted towards the open end (site 



m = 1000) of the sample. It is this latter situation which 
was realized in recent numerical calculations on periodi- 
cally driven hydrogen atoms additionally subjected to an 
intense static electric field [7]. 
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FIG. 3. Survival probability Psurv(t) obtained from distri- 
butions as in Figs. 1 and 2 via Eq. (1). (a) ID liydrogen, initial 
state |no = 100), driving froquoncios and localization param- 
eters uj/2n = 16.45 GHz, £ = 1.0 (full line), 13.16 GHz and 
2.0 (dashed); (b) 3D hydrogen, no = 70, mo = 0, £o = (solid 
line), 15 (dashed), 45 (dotted), lj /2TV = 35.6 GHz, £ = 1.0. 
In this latter case, different -distributions are realized by 
changing the angular momentum io of the initial atomic state, 
always leading to near-algebraic decay (like in ID), though 
with distinct (and non-universal) decay exponents. Only fic- 
titious expansion coefficients Wj = Tj / < F > induce asymp- 
totically "universal" decay Psurv(t) ~ t ' (dash-dotted line), 
7~2 — Q~l.l (indicated by the full line), with Tj, < F >, 



and a ~ 0.9 extracted from the F < 10 
jC = 1.0 data in Fig. i(b). 
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Given the results illustrated in Fig. 2, we can no more 
hope for a universal decay law for the survival probability 
Ps\irv{t) despite the rather universal features of P(r). 
This is finally illustrated in Fig. 3, where we show i-'surv(i) 
for ID and 3D hydrogen in the case of large C: Whilst 
algebraic decay of Psurv(i) ^ is generic (even if mod- 
ulated by local fluctuations [18]), simply due to the broad 
distribution of the wj and Tj as illustrated in Figs. 1 and 
2, the "universal" decay exponent 7 = 2 — a is only ob- 
served in accidental cases - or if we reshuffle the expan- 
sion coefficients of Fig. 2 a) according to Wj rj/(r), 
with (F) the average rate over the entire spectrum, as 
illustrated in Fig. 3 b). In this, however, unrealistic case, 
Ps\iTv{t) asymptotically decays with 7 ~ 1.1, which per- 
fectly matches the prediction 7 = 2 — a, with a = 0.9 
extracted from the small F (< 10"^° a.u.) part of the 
corresponding data set in Fig. 1 b). Therefore, universal 



features of the statistics of the decay rates of classically 
chaotic quantum systems do generically not carry over 
to the time dependent decay of the survival probabil- 
ity of some initial state prepared at an arbitrary loca- 
tion in phase space. This simply expresses the essential 
decorrelatioii of the phase space localization properties of 
the eigenstates of chaotic quantum systems and of their 
asymptotics - which alone determine their decay. 

Support from EU Program QTRANS RTNl-1999- 
08400, is gratefully acknowledged (S.W.). 



[1] 
[2] 

[3] 
[4] 
[5] 

[6; 

[7 

[9: 



[10: 
[11 

[12: 
[13; 

[14 
[15; 

[16: 
[17; 

[is: 



S. Fishman, D. R. Grcmpel, and R. E. Prange, Phys. 
Rev. Lett. 49, 509 (1982). 

G. Casati, I. Guarneri, and D. L. Shepelyansky, IEEE J. 
Quant. Electron. 24, 1420 (1988). 

R.S. MacKay, J.D. Meiss, and I.C. Percival, Phys. Rev. 
Lett. 52, 697 (1984). 

T. Geisel, G. Radons, and J. Rubner, Phys. Rev. Lett. 
57, 2883 (1986). 

M. Reed and B. Simon, Methods of Modem Mathematical 
Physics IV: Analysis of Operators, Academic Press, San 
Diego (1980). 

A. Buchleitner D. Delande, and J.-C. Gay, J. Opt. Soc. 
Am. B 12, 505 (1995). 

G. Benenti et al, Phys. Rev. Lett. 84, 4088 (2000). 
A. Ossipov et al., Phys. Rev. B 64, 224210 (2001). 
A. Buchleitner, Ph.D. thesis, Universite Pierre ct Marie 
Curie, Paris, 1993; A. Buchleitner and D. Delande, 
Chaos, Solitons & Fractals 5, 1125 (1995); Phys. 
Rev. Lett. 75, 1487 (1995); A. Buchleitner et al, 
Phys. Rev. Lett. 75, 3818 (1995). 

R. Ketzmerick et al., Phys. Rev. Lett. 85, 1214 (2000). 
M. Terraneo and I. Guarneri, Eur. Phys. J. B 18, 303 

(2000) . 

P.M. Koch and K.A.H. Leeuwen, Phys. Rep. 255, 289 

(1995) 

A. Krug and A. Buchleitner, Phys. Rev. Lett. 86, 3538 

(2001) . 

G. Casati, G. Maspero, and D. L. Shepelyansky, 
Phys. Rev. Lett. 82, 524 (1999). 

B. Gremaud, D. Delande, and J.-C. Gay, Phys. Rev. Lett. 
70, 1615 (1993). 

G.P. Berman and G.M. Zaslavsky, Phys. Lett. 61A, 295 
(1977); S.W. MacDonald and A.N. Kaufman, Phys. Rev. 
Lett. 42, 1189 (1979); J.G. Leopold and D. Richards, J. 
Phys. B27, 2169 (1994). 

A. Krug, PhD thesis, Ludwig-Maximilians-Universitat 

Miinchen, 2001, available at www.ub.uni- 
muerichen.de / elektronischc_dissertationen / physik/Krug_ 
Andreas.pdf. 

S. Wimberger and A. Buchleitner, J. Phys. A 34, 7181 
(2001). 



4 



